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Abstract:  We  investigate  laser  emission  from  optically-pumped  rota¬ 

tionally  excited  molecular  gases  confined  in  a  metallic  cavity.  To  this  end, 
we  have  developed  a  theoretical  framework  able  to  accurately  describe, 
both  in  the  spatial  and  temporal  domains,  the  molecular  collisional  and 
diffusion  processes  characterizing  the  operation  of  this  class  of  lasers.  The 
effect  on  the  main  lasing  features  of  the  spatial  variation  of  the  electric 
field  intensity  and  the  ohmic  losses  associated  to  each  cavity  mode  are  also 
included  in  our  analysis.  Our  simulations  show  that,  for  the  exemplary  case 
of  methyl  fluoride  gas  confined  in  a  cylindrical  copper  cavity,  the  region 
of  maximum  population  inversion  is  located  near  the  cavity  walls.  Based 
on  this  fact,  our  calculations  show  that  the  lowest  lasing  threshold  intensity 
corresponds  to  the  cavity  mode  that,  while  maximizing  the  spatial  overlap 
between  the  corresponding  population  inversion  and  electric-field  intensity 
distributions,  simultaneously  minimizes  the  absorption  losses  occurring  at 
the  cavity  walls.  The  dependence  of  the  lasing  threshold  intensity  on  both 
the  gas  pressure  and  the  cavity  radius  is  also  analyzed  and  compared  with 
experiment.  We  find  that  as  the  cavity  size  is  varied,  the  interplay  between 
the  overall  gain  of  the  system  and  the  corresponding  ohmic  losses  allows 
for  the  existence  of  an  optimal  cavity  radius  which  minimizes  the  intensity 
threshold  for  a  large  range  of  gas  pressures.  The  theoretical  analysis 
presented  in  this  work  expands  the  current  understanding  of  lasing  action 
in  optically-pumped  far-infrared  lasers  and,  thus,  could  contribute  to  the 
development  of  a  new  class  of  compact  far-infrared  and  terahertz  sources 
able  to  operate  efficiently  at  room  temperature. 
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1.  Introduction 

It  has  been  forty  years  since  the  pioneering  work  of  Chang  and  Bridges  [1]  demonstrated  gen¬ 
eration  of  submillimeter  radiation  via  stimulated  emission  from  optically-pumped  methyl  flu¬ 
oride  (CH3F)  gas.  Since  then,  a  rich  spectrum  of  laser  lines  through  the  50  /im  -  2  mm  wave¬ 
length  region  have  been  obtained  from  a  variety  of  different  molecular  gases  [2-10].  These 
developments  have  made  this  class  of  lasers — usually  referred  to  as  optically-pumped  far  in¬ 
frared  (OPFIR)  lasers — a  central  resource  for  many  scientific  and  technological  applications. 
The  physical  mechanism  underlying  the  operation  of  conventional  OPFIR  lasers  has  also  been 
extensively  studied  in  past  decades  [2-15].  In  essence,  the  gain  in  these  systems  derives  from 
an  optically-driven  population  inversion  between  two  rotational  states  in  an  excited  vibrational 
level  of  a  gas  molecule;  the  laser  emission  frequency  is  determined  by  the  frequency  of  the 
transition  between  these  two  rotational  states.  The  feedback  mechanism  enabling  laser  action 
is  often  provided  by  designing  the  cavity,  in  which  the  molecular  gas  is  confined,  to  act  as  a 
one-dimensional  Fabry-Perot  resonator  for  the  emitted  laser  frequency. 

However,  although  research  on  conventional  OPFIR  lasers  is  generally  considered  as  a  ma¬ 
ture  field,  some  central  challenges  remain  for  the  ultimate  application  of  this  technology,  such 
as  low  efficiency,  poor  tunability,  or  limited  capability  of  integration  into  a  compact  platform. 
More  than  two  decades  ago  it  was  shown  [6]  that  some  of  the  aforementioned  drawbacks  can 
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be  overcome  by  considering  the  active  role  that  high-energy  vibrational  states  play  in  the  laser 
operation  (previous  approaches  assumed  that  the  populations  of  those  vibrational  states  were 
negligible).  These  findings  demonstrated  the  feasibility  of  developing  OPFIR  lasers  able  to  op¬ 
erate  at  large  pressures  (larger  than  the  pressure  cut-off  determined  by  the  so-called  vibrational 
bottleneck  [6,7])  and  featuring  volumes  up  to  four  orders  of  magnitude  smaller  than  those  char¬ 
acterizing  typical  OPFIR  lasers.  The  advent  of  quantum  cascade  lasers  to  replace  CO2  pump 
lasers  may  combine  to  produce  truly  compact  submillimeter-wave  laser  sources  in  the  near 
future. 

Surprisingly,  despite  their  great  promise,  up  to  date  there  has  been  limited  interest  in  further 
exploring  the  potential  of  these  discoveries  [6]  for  the  development  of  a  novel  class  of  com¬ 
pact  far-infrared  and  terahertz  (THz)  sources  able  to  operate  efficiently  at  room  temperature. 
In  order  to  do  that,  it  is  critical  to  develop  a  complete  theoretical  understanding  of  the  physical 
mechanisms  underlying  lasing  action  in  the  considered  systems.  The  purpose  of  this  paper  is 
to  present  such  a  theoretical  analysis.  Specifically,  we  introduce  a  theoretical  framework  which 
generalizes  previous  works  in  this  area  to  allow  for  a  realistic  description  of  the  spatio-temporal 
dynamics  characterizing  the  molecular  collisional  and  diffusion  processes  occurring  in  the  con¬ 
sidered  OPFIR  lasers.  Importantly,  the  effect  on  the  main  laser  features  of  both  the  electric-field 
(ii-field)  intensity  distribution  inside  the  metallic  cell  and  the  ohmic  losses  associated  with  the 
penetration  of  the  fi-field  inside  the  metallic  walls  are  also  included  in  our  theoretical  approach. 
As  we  show  in  this  work,  both  effects  are  particularly  important  for  a  realistic  description  of  the 
operation  of  compact  OPFIR  lasers  based  on  a  metallic  cavity  which  waveguides  the  emitted 
laser  frequency. 

Using  this  theoretical  framework,  we  perform  a  detailed  numerical  analysis  of  laser  emission 
from  optically  pumped  methyl  fluoride  (CH3F)  gas  confined  in  a  cylindrical  copper  cavity,  and 
And  that  the  model  compares  well  against  experimental  results.  Our  calculations  show  that, 
for  this  exemplary  system,  the  region  of  maximum  population  inversion  concentrates  near  the 
cavity  walls.  The  outcome  of  these  calculations  also  implies  that,  as  a  consequence  of  this 
population  inversion  distribution,  the  lowest  threshold  intensity  corresponds  to  the  cavity  mode 
whose  E-fleld  intensity  profile  presents  an  optimal  spatial  overlap  with  the  population  inversion 
distribution,  while  at  the  same  time,  the  ohmic  losses  associated  to  that  cavity  mode  are  min¬ 
imized.  In  the  case  of  the  considered  cylindrical  cavity  both  features  are  satisfied  by  the  TEoi 
mode.  We  also  analyze  how  the  threshold  intensity  depends  on  the  gas  pressure  and  the  cavity 
radius.  The  results  from  these  simulations  are  explained  in  terms  of  the  corresponding  varia¬ 
tion  of  the  overall  gain  of  the  system  and  the  ohmic  losses  associated  with  the  lowest  threshold 
mode.  Remarkably,  we  predict  that  as  the  cavity  size  is  varied,  the  interplay  between  both  fac¬ 
tors  yields  an  optimal  value  of  the  cavity  radius  for  which  the  threshold  intensity  is  minimized 
for  a  large  range  of  gas  pressures. 

The  paper  is  organized  as  follows.  Section  2  discusses  the  theoretical  model  used  throughout 
this  paper.  Section  3  presents  a  detailed  analysis  of  lasing  action  in  the  considered  molecular 
gas  lasers.  Finally,  in  Section  4  we  provide  a  set  of  conclusions  of  this  work. 

2.  Theoretical  framework 

We  start  by  introducing  the  theoretical  model  used  to  describe  the  molecular  dynamics  of  the 
considered  optically-pumped  polyatomic  molecular  gases.  Fig.  1  illustrates  the  energy  level  di¬ 
agram  and  the  different  energy  transfer  mechanisms  that  form  the  basic  structure  of  our  model. 
In  the  following,  we  describe  each  of  the  processes  sketched  in  Fig.  1  and  show  how  to  quantify 
them  using  experimentally  measured  parameters.  Before  starting  this  analysis,  we  point  out  that 
the  theoretical  approach  presented  here  is  an  extension  of  previous  studies  on  methyl  fluoride 
CH3F  [6, 7]  in  order  to  account  for  the  spatial  variations  and  diffusion  of  the  collisional  pro- 
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Fig.  1.  Schematic  diagram  of  the  general  model  used  to  describe  the  dynamics  of  a  OPFIR 
molecular  gas.  Details  on  each  process  labeled  are  described  in  the  text. 


cesses.  We  also  point  out  that  although  for  demonstration  the  model  is  applied  to  the  case  of 
methyl-halide  gases,  our  theoretical  analysis  can  be  easily  generalized  to  most  OPFIR  systems, 
both  pulsed  and  CW,  and  operating  over  a  wide  range  of  experimental  conditions  utilizing  an 
appropriate  gas  that  satisfied  the  desired  pumping  and  lasing  transition  frequencies  [15]. 

The  experimentally-motivated  thermal  pool  concept  has  been  previously  introduced  in 
OPFIR  laser  models  [3,7,  8,  15].  In  the  case  of  molecules,  a  thermal  pool  is  observed  to  be 
a  collection  of  rotational  states  that  are  in  thermal  equilibrium,  and  thus,  are  populated  accord¬ 
ing  to  a  Boltzmann  distribution  at  the  operating  temperature.  Such  equilibrium  may  be  attained, 
for  instance,  due  to  the  fast  relaxation  rates  among  the  levels.  However,  the  total  population  of 
these  pools  varies  through  the  flow  of  energy  governed  by  the  collisional  and  diffusive  pro¬ 
cesses  shown  in  Fig.  1 .  This  fortuitous  consolidation  of  numerous  rotational  levels  into  a  small 
number  of  pools  makes  numerical  simulation  tractable,  and  more  importantly,  makes  rotational 
energy  transfer  studies  possible  in  the  first  place.  In  the  particular  case  of  methyl  halides,  there 
exist  two  thermal  pools  for  each  vibrational  state,  corresponding  to  the  two  symmetry  types 
(A—  and  £— type),  but  without  loss  of  generality  they  will  be  combined  here  to  describe  CW 
OPFIR  laser  operation  at  low  pressure  [8].  In  addition,  all  other  higher  vibrational  states  that 
are  not  connected  by  the  pump  are  also  modeled  as  pools.  This  inclusion  is  important  at  high 
pressure  and  pump  intensity  since  it  provides  an  additional  channel  through  which  molecules 
from  the  excited  vibrational  state  may  be  transferred  out,  and  hence,  reduces  the  absorption  on 
the  corresponding  FIR  lasing  transition  (represented  as  the  pink  arrow  in  the  excited  state  of 
Fig.  1). 

At  the  relatively  low  pressure  regime  considered  in  this  paper,  molecules  are  returned  to  the 
ground  thermal  pool  from  the  pools  of  the  excited  and  higher  excited  states,  shown  as  green 
arrows  in  Fig.  1,  primarily  by  colliding  with  the  walls  of  the  cavity  in  which  the  molecular  gas 
is  confined.  Unlike  molecule-molecule  de-excitation  processes,  these  involve  diffusion  through 
the  gas  towards  the  wall  where  their  rates  of  occurrences  decrease  with  pressure,  and  can  be 
defined  as  kwaii  =  vX^ev/R^  [16].  Amfp  is  the  pressure-dependent  mean  free  path  and  v  is  the 
average  relative  velocity  between  molecules.  At  very  low  pressure  (or  equivalently,  within  Amfp 
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from  the  cavity  walls),  the  transport  mechanism  changes  from  diffusion  to  ballistic  and  the  rate 
becomes  pressure-independent  so  that  A:(iirect  =  t'abs/^j  with  v^bs  dehned  to  be  the  average  ab¬ 
solute  velocity  of  a  molecule.  It  is  clear  from  these  expressions  that  the  wall  rates  are  very 
different  across  the  cavity,  and  we  explicitly  account  for  this  spatial  variation  in  our  numeri¬ 
cal  simulations.  On  the  other  hand,  the  vibrational  relaxation  between  the  excited  and  higher 
excited  pools  can  be  dominated  by  molecule-molecule  collisions,  in  which  case,  the  transition 
rates  between  them  are  ^exc  =  Mott'<7exc-  Here,  A^tot  is  the  total  density  of  molecules  and  ffexc  is 
the  collision  cross  section  measured  from  experiments. 

Non-equilibrium  among  rotational  and  vibrational  states  in  the  model  is  solely  introduced  by 
the  pump  (i.e.  the  blue  arrow  in  Fig.  1),  and  the  resulting  nonthermal  distribution  of  population 
in  some  of  the  rotational  levels  can  no  longer  be  considered  part  of  the  thermal  pool  but  must 
be  modeled  individually.  These  nonthermal  states,  illustrated  in  Fig.  1  as  level  Ji,  72-  ■  ■  etc.  in 
each  vibrational  state,  only  exist  in  states  that  are  of  the  same  symmetry  type  and  vibrational 
state  as  those  connected  by  the  pump.  The  total  population  of  any  state  7,  (i.e.  Mot.y,)  is  simply 
the  sum  of  its  nonthermal  and  thermal  parts,  Mot.y,  =  Nj-  +  /y,Npooi,  where  fj-  is  the  fraction 
of  the  total  pool  population  Npooi  in  state  7,  and  Nj-  is  the  amount  the  state  population  differs 
from  rotational  thermal  equilibrium  because  of  the  pump.  Note  that  Nj.  can  be  positive  (excited 
vibrational  state)  or  negative  (ground  vibrational  state). 

In  general,  the  pump  will  create  signihcant  non-equilibrium  not  only  in  the  states  that  it 
directly  connects,  but  also  in  those  adjacent  to  it.  To  account  for  this,  our  numerical  model  al¬ 
lows  one  to  include  as  many  rotational  states  as  are  necessary  to  model  such  a  pump-induced 
non-equilibrium  accurately.  Generally,  the  number  of  J  states  to  include  depends  on  the  ratio 
of  the  dipole-dipole  rate  to  the  K-swap  rate  (see  end  of  paragraph  for  the  description  of  these 
two  processes)  —  the  slower  the  K-swap  rate,  the  more  J  levels  one  needs  to  consider.  As  al¬ 
ready  pointed  out,  the  total  population  of  any  nonthermal  state  in  this  model  is  the  combined 
population  from  both  the  nonthermal  and  thermal  contributions,  but  since  the  thermal  pools 
are  in  rotational  equilibrium,  the  rotational  relaxations,  shown  as  black  arrows  in  Fig.  1,  are 
only  modeled  for  the  nonthermal  population.  In  methyl  fluoride,  the  rotational  nonthermal  pro¬ 
cesses  are  primarily  dipole-dipole  collisions  whose  rates  are  computed  from  the  dipole-dipole 
collision  cross  section  Odd  measured  from  experiments  [9].  Note  that  their  allowed  transitions 
are  restricted  to  zi/  =  ±  1  at  the  same  K  =  Kj  as  the  pumping  transition.  Lastly,  the  K-swap 
process  [8]  (i.e.  the  red  arrows  in  Fig.  1)  allows  population  to  be  transferred  out  of  nonthermal 
states  into  the  rest  of  the  rotational  levels  within  the  same  thermal  pool,  and  its  rate  constant 
can  again  be  measured  experimentally.  The  K-swap  collisional  processes  (which  encompass  the 
AK=3n  process  and  the  V-swap  process  [8, 15]  when  the  A—  and  £— symmetry  thermal  pools 
are  combined  as  they  are  here)  convert  the  nonthermal  molecules  to  thermal  molecules  but  not 
vice  versa,  and  thus,  allowing  the  nonthermal  states  to  return  to  zero  when  the  pump  is  turned 
off. 

All  of  the  aforementioned  energy  transfer  mechanisms  except  K-swap  include  reverse  pro¬ 
cesses  which  account  for  the  appropriate  degeneracy  of  the  upper  and  lower  levels  so  that 
microscopic  reversibility  is  satished.  Moreover,  spatial  diffusion  of  the  molecules  within  the 
cavity,  though  a  relatively  slow  process,  is  also  included  in  order  to  account  for  the  molecular 
gas  dynamics  within  the  cavity.  The  pump  transition  rate,  assuming  a  cylindrical  waveguide 
cavity,  is  determined  from  R^ump  =  <xPpump/^R^Ntotfiov/hViR,  where  Fpump  is  the  pump  power, 
/low  =  fji  is  the  fraction  of  molecules  in  the  lower  ro-vibrational  state  7i  (as  depicted  in  Fig.  1), 
ViR  is  the  pump  infrared  frequency,  and  a,  which  takes  into  account  both  Doppler  broadening 
and  mismatch  between  pump  and  absorption  lines,  measures  the  gas’  ability  to  absorb  the  in¬ 
frared  radiation.  Hence,  pump  transition  rate  may  be  increased  by  reducing  R,  increasing  Fpump, 
or  decreasing  Vir.  The  last  of  these,  however,  is  constrained  by  the  need  to  match  the  frequency 
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of  the  laser  line  with  that  of  a  ro-vibrational  transition  in  the  molecules,  and  in  general,  cannot 
be  varied  at  will.  Using  such  an  implementation,  our  model  accommodates  both  CW  and  pulsed 
pumping  with  arbitrary  spatial  profiles. 

Lastly,  we  describe  the  various  processes  in  the  model  via  rate  equations,  within  the  rate- 
equation  approximation  [17].  The  most  general  rate  description  of  a  rotational  level  would  in¬ 
clude  rotational  relaxation  terms  due  to  dipole-dipole  collisions  among  the  non-thermal  J  levels, 
terms  representing  the  pumping  mechanism,  K-swapping  terms  that  allow  flow  of  molecules  out 
of  the  non-thermal  levels  (with  equilibrating  effects  on  non-thermal  levels),  and  the  stimulated 
emission  terms  at  the  FIR  wavelength.  Since  we  are  interested  in  the  steady  state  distribution  of 
molecules  across  the  cavity,  spatial  diffusion  terms  are  also  added  to  the  rate  equations  describ¬ 
ing  each  levels.  Stimulated  and  spontaneous  transitions  at  the  lasing  wavelength  can,  however, 
be  neglected  in  the  rate  equations  when  operated  near  the  threshold  if  the  spontaneous  emission 
lifetime,  Tspont,  is  relatively  longer  than  all  other  transition  lifetimes  in  the  considered  system.  In 
methyl  fluoride,  for  instance,  Tspont  is  on  the  order  of  lOi  or  more.  In  contrast,  the  next  longest 
lifetime  in  the  model  occurs  for  the  vibrational  relaxation  terms  and  yet,  has  lifetime  on  the 
order  of  1  x  lO^^i  or  less.  The  relative  magnitude  of  these  radiative  terms  involving  Tspont  in 
the  rate  equations  are  further  diminished  by  the  fact  that  only  a  fraction  of  these  transitions 
contributes  to  the  desired  laser  mode,  especially  when  the  cavity  is  large.  With  that,  the  rate 
equation  describing  the  molecular  density  of  J2  level  in  the  excited  vibrational  state  (i.e. 
of  Fig.  1  is  given  by 


Jt 


3  ,  2  ^exc/j. 

--Rpump(r,f)AA?pump(r,f)  X  —  +  S  ^ 
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where  AA^pump(r,f)  =  A^toU2(*’>0  with  being  the  total  molecular 

density  of  the  level  Ji  in  the  ground  vibrational  state  after  accounting  for  the  nonthermal  density 
and  the  contribution  from  the  corresponding  fraction  of  the  thermal  pool,  and  =  27  -f  1  is 
the  degeneracy  of  the  ground  state  level.  Similar  definitions  hold  for  and  g}-.  Tj,,i 

for  !  =  1,2,3  is  the  decay  lifetime  out  of  to  the  two  adjacent  rotational  levels  and  also 
the  thermal  pool,  for  j  =  1,2  represents  the  decay  into  from  its  two  adjacent  levels 
and  Finally,  D  is  the  diffusion  constant  given  by  the  product  of  the  molecule  velocity 
and  mean  free  path.  In  the  same  way,  the  thermal  pool  equation  has  terms  from  the  K-swap 
and  spatial  diffusion  processes,  as  well  as  vibrational  relaxations  due  to  molecule-molecule 
collisions  and  diffusion  to  wall.  The  equation  describing  the  excited  vibrational  pool  in  Fig.  1 
is  then  given  by 


(r,0 


dt 


(2) 


The  first  two  terms  relate  to  the  vibrational  relaxations  among  the  thermal  pools,  where  Texc,/t(r) 
for  ^  =  1 , 2  is  the  decay  lifetime  from  the  excited  thermal  pool  to  the  ground  and  higher  excited 
pools,  while  T„^exc(r)  is  the  reverse.  As  noted  earlier,  the  values  of  these  diffusive  lifetimes 
depend  on  their  distances  from  the  cell  wall.  T™^^  for  m=  1,2, .. .  ,mo  in  the  third  term  is  the 
lifetime  of  the  transition  from  each  of  the  mo  rotational  levels  to  the  excited  pool  due  to  the 
K-swap  molecule-molecule  collision  process.  As  in  Eq.  (1),  the  last  term  represents  spatial  di¬ 
ffusion.  Similar  equations  to  Eq.  (1)  and  Eq.  (2)  are  used  to  represent  every  nonthermal  level  or 
thermal  pool  that  contributes  to  the  lasing  action,  which  in  turn  depends  on  their  degree  of  non¬ 
equilibrium  caused  by  the  pump  (see  Appendix).  This  way,  the  complete  set  of  rate  equations 
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describing  our  laser  system  obeys  the  diffusion  equation  dNioi{r,t)/dt  =  DV^Niot(j,t),  and  so, 
conserves  the  total  molecular  density. 

In  our  simulations,  we  consider  only  the  immediate  nonthermal  levels  adjacent  to  the  pumped 
levels  and  initialize  the  density  of  molecules  in  the  thermal  pools  according  to  a  Boltzmann- 
distributed  partition  based  on  their  respective  vibrational  energies.  The  non-thermal  7,  levels  are 
set  equal  to  zero  before  any  pumping  occurs.  The  set  of  equations  is  then  evolved  in  time  until 
steady  state  is  reached,  allowing  the  temporal  development  of  the  lasing  action  to  be  tracked. 

Laser  action  in  this  class  of  OPFIR  molecular  gas  lasers  is  obtained  as  follows:  Pumping 
is  achieved  with  a  line-tunable  laser  (e.g.  CO2  laser)  which  excites  molecules  from  a  specific 
rotational  level  in  the  ground  vibrational  state  into  a  specific  rotational  level  in  an  excited  vibra¬ 
tional  state.  Because  the  excited  vibrational  level  is  comparatively  empty,  these  photo-excited 
molecules  create  a  population  inversion  between  the  pumped  rotational  state  and  the  one  im¬ 
mediately  below  it.  This  leads  directly  to  stimulated  emission  between  rotational  states  in  the 
excited  vibrational  state.  These  photo-excited  non-equilibrium  molecules  are  subsequently  ro- 
tationally  and  vibrationally  relaxed  to  the  ground  state  through  collisions  with  other  molecules 
and  with  the  chamber  walls,  respectively. 

The  population  inversion  responsible  for  lasing  exists  as  long  as  the  pump-induced  difference 
in  nonthermal  population  between  the  upper  and  lower  state  Nj.  —Nj--i  is  greater  than  the  dif¬ 
ference  in  the  corresponding  pool-contributed  thermal  populations  (/y.  — /y,-i)7fpooi  for  those 
states.  As  the  pressure  grows,  the  pool  molecules  have  an  increasingly  difficult  time  reaching 
the  walls  and  de-exciting;  thus,  the  thermal  contribution  from  the  pool  increasingly  quenches 
the  non-thermal  population  inversion.  Pool  quenching  creates  the  vibrational  bottleneck  that 
was  once  believed  to  limit  OPFIR  laser  operation  to  low  pressures  regardless  of  pump  inten¬ 
sity.  However,  it  is  now  understood  that  as  the  excited  vibrational  state  thermal  pool  fills,  near 
resonant  collisions  among  its  constituents  can  excite  molecules  to  even  higher-lying  vibrational 
levels,  thereby  providing  an  alternate  relaxation  pathway  [6].  Therefore,  the  pressure  at  which 
vibrational  bottleneck  occurs  for  a  given  laser  geometry  depends  on  the  pump  intensity.  By 
concentrating  on  low  pressure  operation  in  this  manuscript,  these  additional  pathways  may  be 
ignored  to  simplify  the  analysis  without  affecting  the  generality  of  the  results. 

When  molecules  are  photo-excited  from  the  ground  vibrational  state  to  an  excited  state, 
two  types  of  population  inversion  can  be  induced,  and  lasing  arises  from  each.  The  first  is  the 
primary  lasing  transition  already  described  between  the  pumped  rotational  state  in  the  excited 
vibrational  level  and  the  rotational  state  below  it.  The  second  is  the  refilling  transition  between 
the  IR  pump  depleted  rotational  state  in  the  ground  vibrational  level  and  the  rotational  state 
above  it.  When  the  pump  is  eventually  turned  off,  the  K-swap  process  allows  equilibrium  con¬ 
ditions  in  the  molecules  to  be  restored.  In  the  following  illustration  of  this  theoretical  approach, 
without  loss  of  generality,  we  present  results  only  for  the  primary  lasing  transition. 

3.  Results  and  discussion 

Figure  2(a)  illustrates  the  setup  of  the  OPFIR  laser  system  considered  in  this  section.  A  cylin¬ 
drical  waveguide  resonator  of  length  Lceii  and  radius  R  is  filled  with  a  suitable  gas  that  lases  at 
the  desired  THz  frequency,  while  pump  power  at  a  much  higher  IR  frequency  enters  the  system 
from  the  front  window.  Depending  on  the  reflectivities  of  the  front  and  back  windows  at  THz, 
lasing  output  power  can  then  escape  from  the  cavity  via  both  channels.  In  order  to  compare 
our  numerical  calculations  with  past  experimental  results,  the  cavity  modeled  in  the  follow¬ 
ing  simulation  has  length  Lceii  =  12  cm  and  front  window  reflectivity  Fftont  =  0.96.  The  back 
window  is  opaque  (i.e.  Fback  =  !)■  Furthermore,  R  is  varied  across  a  range  of  commercially 
available  values  from  the  smallest  possible  before  modal  cutoff  at  0.08  cm  to  1  cm.  The  cavity 
is  then  filled  with  an  isotopic  isomer  of  methyl  fluoride  gas,  ^^CH3F,  where  a  total  of  six  rota- 
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Fig.  2.  (a)  The  OPFIR  laser  system  considered  in  our  numerical  modeling:  A  cylindrical 
waveguide  resonator  of  length  Lceii  and  radius  R  is  filled  with  a  suitable  gas  that  lases  at 
the  desired  TFlz  frequency,  while  pump  power  at  a  much  higher  IR  frequency  enters  the 
system  from  the  front  window.  In  our  example,  *^CH3F  is  pumped  with  CO2  laser  at  31 
THz  to  produce  lasing  at  0.245  THz.  (h)  Frequency  dependence  of  the  ohmic  losses  for  the 
five  lowest  order  modes  (i.e.  TEn,  TMqi,  TE21,  TEqi,  TMh)  and  three  higher  order  ones 
with  low  losses  (i.e.  TE02,  TE12,  TE22).  Plot  is  for  a  copper  cavity  with  R  =  0.26  cm  (and 
assuming  Leell  ^  so  that  several  modes  (not  all  shown  in  plot)  have  cut-offs  helow  the 
THz  lasing  frequency.  Despite  this,  only  the  lowest  loss  mode  TEqi  exist  in  the  cavity  near 
the  lasing  threshold.  Vertical  green  line  shows  the  hand  of  frequency  at  which  our  system 
operates,  (c)  Intensity  profile  of  the  three  lowest  loss  modes  supported  by  the  R  =  0.26  cm 
waveguide  cavity. 


tional  levels,  7vo  =3,4,5  in  the  ground  vibrational  state  and  =4,5,6  in  the  excited  state, 
are  included  (see  detailed  rate  equations  in  the  Appendix).  Here,  the  K-swap  processes  occur  at 
roughly  1-2  times  the  rates  of  the  rotational  transitions,  and  through  this  process,  the  lower  and 
higher  adjacent  J  levels  can  be  modeled  as  part  of  their  respective  pool.  For  low  pressure  op¬ 
eration,  we  include  the  doubly  degenerate  Ve  state  as  our  only  higher  excited  vibrational  state. 
As  we  shall  see  next,  these  relatively  small  number  of  rotational-vibrational  levels  are  already 
sufficient  to  attain  reasonable  match  between  numerical  predictions  and  experimental  data.  The 
system  is  uniformly  pumped  at  room  temperature  {T  =  300  K)  with  CW  CO2  laser  at  31  THz 
in  which  molecules  are  excited  between  7vo  =  4  and  Jy^  =  5  to  produce  primary  lasing  tran¬ 
sition  at  0.245  THz  in  Tys  =4  —  5  of  the  V3  vibrational  state,  with  Tspom  roughly  chosen  to  be 
15i  [19].  It  is  also  noted  that  the  occupancy  of  V3  is  comparatively  low  (/  =  k.  1/150), 

making  the  condition  favorable  for  the  creation  of  population  inversion  between  Jy^  =  5  and  the 


#142492 -$15.00  USD 
(C)  2011  OS  A 


Received  10  Feb  2011;  accepted  21  Mar  2011;  published  4  Apr  201 1 
1 1  April  2011/  Vol.  19,  No.  8  /  OPTICS  EXPRESS  7520 


level  immediately  below.  Lastly,  Ki  in  Fig.  1  equals  3  for  *^CH3F  and  we  have  approximated 
the  K-swap  process  as  the  faster  A/ir=3n  process  where  the  V-swap  rate  is  slow  and  primarily 
equilibrates  A  and  E  states  in  the  pool.  Other  molecules  with  a  slower  or  absent  A/ir=3n  process 
will  need  to  approximate  the  K-swap  process  by  the  V-swap  process.  The  rate  constants  and 
other  salient  molecular  parameters  used  are  provided  in  the  Appendix. 

Figure  2(b)  shows  a  plot  of  the  ohmic  losses  in  a  /?  =  0.26  cm  copper  waveguide  cavity  as  a 
function  of  frequency.  This  loss  prediction  is  made  assuming  Lceii  ^  Rso  that  only  the  circular 
plane  of  the  cavity  matters  and  that  the  copper  conductivity  is  high  enough  that  the  wall  currents 
flow  uniformly  within  a  skin  depth  of  the  surface  [18].  At  the  operating  frequency,  shown  as  a 
green  vertical  band  in  figure,  the  cavity  supports  multiple  modes  but  only  the  five  lowest  order 
modes  and  three  higher  order  ones  that  have  low  ohmic  losses  are  included  in  the  plot.  The 
cut-off  frequency  of  each  of  these  modes  can  also  be  deduced  from  the  figure.  The  lowest  loss 
mode,  TEoi,  will  be  the  first  to  lase  upon  threshold  and  remains  the  only  mode  present  for 
pump  power  near  the  threshold,  despite  the  cavity’s  ability  to  support  multiple  modes  at  THz. 
Alternatively,  one  can  ensure  single-mode  operation  even  for  large  pump  power  by  selecting  the 
smallest  possible  cavity  size  with  radius  ^  /Ithz/2,  where  Athz  is  the  wavelength  of  the  THz 
lasing  output.  For  our  purpose,  it  suffices  to  simply  operate  near  the  threshold  for  single-mode 
operation  since  this  also  facilitates  independent  design  optimizations  of  our  system  based  on 
molecular  gas  physics  and  photonics  considerations.  Figure  2(c)  depicts  the  intensity  profiles 
of  the  three  lowest  loss  modes  of  the  cavity.  It  may  be  noted  that  the  magnitude  of  ohmic  loss 
scales  roughly  with  the  fraction  of  modal  intensity  residing  near  the  cell  wall,  which  in  turn 
relates  to  the  amount  of  penetration,  to  within  a  skin  depth,  into  the  wall.  At  low  enough  pump 
rate,  only  losses  of  TEqi  mode  (both  ohmic  losses  and  leakage  through  the  end  mirrors)  can  be 
compensated  by  the  gain  produced  in  *^CH3E  while  the  other  modes  remain  suppressed.  Thus, 
we  shall  consider  TEqi  to  be  the  lowest-threshold  lasing  mode. 

Eigure  3  plots  typical  spatio-temporal  results  of  the  numerical  simulation.  In  particular,  the 
evolution  of  the  laser  transition  AN (Athz)  and  the  V3  thermal  pool  to  steady  state  are  presented 
and  compared  at  250  and  350  mTorr,  assuming  Ppump  to  be  10  watts  in  CW  operation.  A  few 
initial  observations  are  in  order:  (i)  the  radial  spatial  profile  of  AN{r)  is  the  inverse  of  that  for 
V3  pool.  This  is  the  pool  quenching  that  leads  to  the  vibrational  bottleneck  mentioned  earlier. 
Note  the  rapid  removal  of  molecules  from  V3  near  the  cell  wall  is  critical  to  maintaining  the 
inversion,  (ii)  The  cell  center  takes  a  long  time  to  attain  steady  state  given  that  the  wall  rate 
is  slow  far  away  from  the  wall  (i.e.  molecules  have  to  travel  a  long  distance  to  reach  the  cell 
boundary).  (Hi)  Similarly,  a  longer  time  is  required  to  reach  steady  state  at  a  higher  pressure  due 
to  the  decreased  near  the  cell  center.  In  Eig.  3,  steady  state  is  attained  after  100  at  250 
mTorr  while  it  takes  300  jJLs  to  do  so  at  350  mTorr.  ( iv)  The  magnitude  of  the  variation  across  the 
radial  direction  of  the  cavity  is  greater  at  higher  pressure,  as  will  be  explained  shortly.  While  we 
mostly  deal  with  CW  systems  in  this  paper,  the  ability  to  track  the  full  temporal  development 
of  the  laser  action  allows  pulsed  systems  to  be  studied  as  well. 

We  next  verify  our  numerical  model  against  experimental  results  for  a  /?  =  0.26  cm  copper 
cavity  ^^CH3E  laser.  The  experimental  configuration  was  identical  to  that  described  in  [6], 
with  a  1  mm  diameter  output  coupler  hole  in  the  front  reflector  through  which  the  IR  pump 
and  THz  laser  power  were  transmitted  [Eig.  2(a)].  The  slope  efficiency  was  measured  as  the 
ratio  of  the  THz  laser  power  emitted  to  the  corresponding  IR  pump  intensity,  and  the  threshold 
pump  intensity  was  the  extrapolation  of  this  slope  to  zero  THz  power.  Numerical  predictions 
and  experimental  measurements  of  both  the  pump  power  and  slope  efficiency  at  threshold  are 
compared  at  two  pressures:  100  mTorr  and  300  mTorr.  Here,  we  adopt  a  semi-analytic  approach 
for  these  comparisons,  where  knowledge  of  the  losses  in  the  OPEIR  system  and  population 
inversions  calculated  with  our  model  are  used  simultaneously.  In  the  following,  we  first  discuss 
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Fig.  3.  Time  dependent  properties  of  laser  system  with  CW  pumping  at  /pump  =  10  W.  (a) 
Time  evolution  of  the  inversion  in  V3  state  at  250  mTorr  across  the  radial  axis  of  the  cavity 
until  steady  state  behavior  is  observed,  (b)  Time  evolution  of  the  V3  thermal  pool  population 
at  250  mTorr  across  the  radial  axis  of  the  cavity  until  steady  state  behavior  is  observed,  (c) 
Same  as  in  (a)  except  operated  at  350  mTorr  (using  the  same  magnitude  range  for  the  color 
bar),  (d)  Same  as  in  (b)  except  operated  at  350  mTorr  (using  the  same  magnitude  range  for 
the  color  bar). 


the  match  for  the  threshold  power  intensity  before  addressing  that  of  the  slope  efficiency. 

The  blue  lines  in  Figs.  4(a)  and  4(b)  depict  the  unsaturated  effective  gain  values  of  the  system, 
as  predicted  by  our  model,  for  a  range  of  pump  intensity  at  100  and  300  mTorr  respectively. 
The  gain  for  the  7  =  4  —  5  lasing  transition  in  V3  state  of  *^CH3F  may  be  written  as  [20, 21] 


(7(Athz))  =  (AAi) 


■^Hz 

4  Tspont  A  T  (  AtHz  ) 


(3) 


where  (AN)  is  the  pressure-dependent  steady-state  effective  inversion  defined  as 

|2  , 


{AN{Q)  = 


(4) 


Eo(r)  is  the  normalized  mode  profile  of  the  passive  cavity  (/  dreon^(r)  |Eo(r)p  =  1)  and  Fact 
being  the  volume  of  active  region  of  the  considered  structure.  In  our  case.  Fact  is  the  entire 
cavity  (note  that  this  implies  that  the  gas’s  behavior  is  less  relevant  in  regions  where  the  modal 
intensity  is  small  but  plays  an  important  role  in  regions  of  high  field  intensity).  In  the  considered 
case,  AN (r)  =  (r)  —  j=4i^) /s^Ls  with  N^^^J  being  the  total  molecular  density 

of  the  upper  (7  =  5)  and  lower  (7  =  4)  lasing  transition  level  in  V3  state  after  accounting  for 
the  fractional  contribution  from  the  excited  thermal  pool,  and  gp  =  27  -F  1  is  the  degeneracy 
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of  the  excited  state’s  J‘^  level.  Eo(r)  is  the  normalized  field  of  the  TEqi  mode  and  n  Ri  1  for 
the  gas.  The  width  of  the  transition  line,  Av,  is  broadened  due  to  a  combination  of  Doppler 
effects  and  molecular  collisions  and  has  a  strong  dependence  on  both  pressure  and  transition 
frequency.  The  former  is  associated  with  the  Gaussian  lineshape  while  the  latter,  which  leads 
to  pressure  broadening,  is  Lorentzian.  Hence,  Av  =  (Av^  +  AVp)^/^  where  the  two  terms  are 
the  Doppler  and  pressure  broadening  widths  respectively,  and  the  overall  lineshape  function 
used  is  the  one  that  corresponds  to  the  larger  broadening  parameter  of  the  two.  It  is  noted 
that  for  the  same  half  width,  pressure-broadened  lines  are  stronger  far  away  from  the  line- 
center  than  the  Doppler-broadened  ones.  The  gain  defined  in  Eq.  (3)  accounts  for  the  number 
of  modes  only  in  the  frequency  interval  Vthz±Av/2  by  assuming  a  large  optical  resonator 
whose  dimension  is  several  wavelengths  so  that  its  exact  shape  becomes  less  significant.  We 
find  that  this  applies  in  our  case  since  the  cylindrical  copper  waveguide  is  almost  IOOAthz  long 
and  a  few  AthzS  in  diameter.  Moreover,  the  inversion  density,  and  hence  gain,  is  not  saturated 
by  stimulated  emission  for  the  entire  range  of  input  intensity  plotted.  As  we  are  only  interested 
in  the  threshold  value  in  this  part  of  the  work,  such  approach  remains  valid. 

The  threshold  for  lasing  may  be  approximated  (neglecting  spontaneous  emission)  to  occur 
when  the  gain  balances  the  total  losses  of  the  OPEIR  system.  Losses  in  our  setup  are  assumed 
to  be  due  to  either  ohmic  losses  from  the  E-field  penetration  into  the  copper  wall  or  to  energy 
leakage  from  the  front  mirror.  Eor  Lceii  ^  R,  we  can  treat  the  system  as  one-dimensional  along 
the  cylinder  axis  for  loss  analysis  and  quantify  its  coefficient  as  2aohmic  —  In(rfrontrback) /LceU, 
where  aohmic  is  the  amplitude  loss  coefficient  shown  in  Eig.  2(b).  The  magnitudes  of  these 
loss  coefficients  are  indicated  as  horizontal  red  lines  in  Eigs.  4(a)  and  4(b).  As  such,  the  input 
threshold  intensity  predicted  from  our  model  is  determined  at  the  intersection  between  the 
red  and  blue  lines,  which  indicates  the  point  at  which  losses  are  compensated  by  the  gain  at 
■^THz  so  that  any  further  increase  in  pump  intensity  saturates  the  gain  and  is  thereby  channeled 
into  coherent  output  power  that  escapes  from  the  front  window.  Erom  Eigs.  4(a)  and  4(b),  the 
numerical  thresholds  attained  are  0.78  and  2.42  W/cm^  for  the  laser  operated  at  100  mTorr  and 
300  mTorr,  respectively. 

To  make  realistic  comparisons  with  experiment,  the  decreasing  absorption  of  the  IR  pump 
beam  as  it  propagates  along  the  cavity  has  to  be  considered.  The  pressure-dependent  IR  ab¬ 
sorption  coefficient,  a  ,  for  *^CH3E  is  71  m^'Torr^*  [14]  so  that  the  ability  of  the  pump  to 
excite  the  gas  decreases  with  increasing  pressure.  The  effect  is  to  reduce  the  amount  of  IR 
power  effectively  pumping  the  gas  and  can  be  captured  by  defining  an  incident  pump  power 
^pump  =  Tpump/K",  where  tc  is  a  pressure-dependent  factor  that  accounts  for  the  net  fraction  of 
incident  IR  power  available  for  excitation  in  the  Lceii  long  cavity.  This  factor  is  simply  given  by 

(1/Lceii)  X  ^dx.  Hence,  K  =  0.67  and  0.36  for  our  100  mTorr  and  300  mTorr  setups 

resulting  in  actual  pump  thresholds  (Epu^p/;rR^)  of  1.2  and  6.7  W/cm^,  respectively.  These  re¬ 
sults  compare  well  with  experimentally  measured  thresholds  of  1 . 1  ±  0.4  and  5 .5  ±  1 .5  W/cm^. 
Possible  sources  of  discrepancy  in  the  comparison  include  instabilities  that  make  measuring 
the  threshold  difficult,  deviation  from  spatially  uniform  pumping  in  the  experiment,  the  one¬ 
dimensional  treatment  in  balancing  gain-loss,  and  the  possibility  that  the  emission  was  from 
the  refilling  transition  in  the  Vq  state.  Even  so,  the  current  treatment  provides  adequately  ac¬ 
curate  predictions  of  the  lasing  thresholds.  More  importantly,  this  comparison  acts  to  verify 
our  claim  that  the  aforementioned  model  of  the  OPEIR  is  physically  robust  and  rigorous,  yet 
numerically  tractable. 

Next,  we  present  a  similar  comparison  for  the  slope  efficiency  at  threshold  using  a  semi- 
analytic  approach.  The  set  of  rate  equations  in  Eq.  (1),  with  the  inclusion  of  the  induced  lasing 
transition  rate,  are  solved  analytically  at  steady-state  to  obtain  the  following  approximate  rela- 
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Fig.  4.  (a)  Method  of  threshold  predictions  via  gain-loss  balancing  used  in  plotting  (c)  for 
every  pressure  point.  Flere,  we  operate  at  100  mTorr  with  R  =  0.26  cm.  The  blue  line  is  the 
unsaturated  effective  gain  value  for  a  range  of  pump  intensity  predicted  by  the  aforemen¬ 
tioned  model  while  the  red  line  is  the  magnitude  of  the  loss  coefficient  due  to  absorption 
by  the  metallic  cavity  and  energy  leakage  from  the  front  window.  Intersection  between  the 
two  lines  corresponds  to  the  lasing  threshold,  (b)  Same  as  in  (a)  except  that  the  pressure 
is  at  300  mTorr.  (c)  Left  axis  depicts  the  relationship  between  the  threshold  intensity  and 
pressure  for  a  range  of  R  values  from  0.08  cm  to  1  cm  as  predicted  from  numerical  model. 
The  corresponding  linewidth,  Av,  of  the  lasing  transition  (over  which  it  may  be  tuned)  is 
also  illustrated  in  the  top  axis.  A  general  trend  exists  such  that  the  lasing  threshold  increases 
with  pressure.  The  right  axis  shows  the  mean  free  path,  Amfp.  within  the  gaseous  system 
which  at  a  fixed  temperature,  is  inversely  proportional  to  pressure.  Inset  is  the  ohmic  loss 
[see  Fig.  2(b)]  of  the  TEqi  mode  at  0.245  TFlz,  with  the  radii  of  interest  marked  as  square 
markers. 


tionship  between  output  power  ^out  and  input  power  Pin  above  threshold; 


Pout  =  /outAat  (  ^  -  1  )  (5) 

where  /out  accounts  for  the  fraction  of  emitted  power  that  escapes  via  the  front  window,  and 
is  the  IR  power  at  threshold.  Psat  is  the  saturation  power  at  THz  in  which  the  gain  in  the  system 
is  half  its  unsaturated  value.  In  this  particular  system,  Psat  =  JiO>iHz  /v^ct  / Teff  where 

AA^th  is  the  inversion  at  threshold  and  Tgff  is  an  effective  lifetime  that  accounts  for  the  various 
rates  between  the  rotational  levels  in  the  lasing  V3  state,  and  is  found  to  be 
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The  lifetimes  in  Eq.  (6)  describe  the  transition  rates  in  the  excited  V3  state  where  lasing  occurs 
and  are  defined  in  accordance  to  their  definitions  in  Fig.  1  and  Eq.  (1).  With  that,  the  slope 
dPom.ldPm  is  /outEsat/ffn  ■  Using  numerical  results  of  AA^th  and  fj*  found  from  the  full  model  and 
accounting  for  the  pressure-dependent  IR  absorption  discussed  in  the  previous  paragraph,  we 
obtain  a  threshold  slope  efficiency  of  24.5  x  10^®  and  39.2  x  10^®  at  100  mTorr  and  300  mTorr, 
respectively.  Note  that  the  best  slope  efficiency  for  the  considered  *^CH3F  laser  is  8  x  10^^, 
which  is  the  Manley-Rowe  limit  when  the  quantum  efficiency  is  1.  These  calculated  values 
compare  reasonably  well  with  experimental  slope  values  of  (44  ±22)  x  10^®  and  (48  ±24)  x 
10^®,  for  which  the  greatest  source  of  uncertainty  is  the  calibration  of  the  the  Golay  cell  used 
to  detect  the  THz  radiation.  Given  this,  these  results  further  validate  our  model. 

We  extend  the  threshold  determination  method  of  Figs.  4(a)  and  4(b)  to  study  how  the  thresh¬ 
old  intensity  behaves  as  both  cell  radius  and  pressure  are  varied  in  Fig.  4(c).  A  general  trend  ex¬ 
ists  such  that  the  lasing  threshold  increases  with  pressure  for  every  radius  plotted.  This  can  pri¬ 
marily  be  attributed  to  the  fact  that  the  linewidth,  Av(/Ithz),  broadens  with  increasing  pressure 
so  that  (7(Athz))  reduces  [see  Eq.  (3)].  Note  that  the  variation  of  the  linewidth  with  pressure  is 
captured  in  the  top  axis  of  the  figure.  Since  the  total  loss  of  the  system  is  pressure-independent, 
such  gain  reduction  increases  the  threshold  linearly  with  pressure  for  each  value  of  R.  Second 
order  effects  that  further  raise  the  threshold  at  high  pressure  originate  from  the  vibrational  bot¬ 
tleneck:  A  higher  pressure  impedes  the  molecular  diffusion  to  the  cell  wall  and  the  subsequent 
decrease  in  wall  collision  rates  reduces  the  flow  of  molecules  out  of  the  the  excited  thermal 
pool,  leading  to  enhanced  absorption  and  hence,  a  weaker  inversion  (AN)  and  gain.  Equiv¬ 
alently,  the  same  effects  on  (AN)  can  be  understood  by  recognizing  that  Amfp  progressively 
makes  up  a  smaller  region  near  the  cell  wall  as  the  pressure  increases  such  that  a  greater  pro¬ 
portion  of  the  molecules  are  in  the  non-ballistic  regime,  which  also  leads  to  slower  kwaii-  Red 
dashed  line  plotted  against  the  right  axis  depicts  the  pressure  dependence  of  Amfp-  Over  the 
range  of  low  pressures  considered,  the  vibrational  bottleneck  is  only  observed  for  large  radius 
cells  (i.e.  R  >  0.26  cm  in  the  current  setup)  while  for  smaller  cavities,  there  is  no  vibrational 
bottleneck  and  the  threshold  dependence  on  pressure  remains  approximately  linear. 

From  Fig.  4(c),  it  can  also  be  seen  that  the  threshold  intensity  increases  with  radius  for  R 
values  larger  than  0.26  cm  but  decreases  with  radius  for  R  values  smaller  than  0.26  cm.  To 
understand  why  R  —  0.26  cm  is  nearly  optimal,  the  reader  is  referred  to  the  inset  that  illustrates 
the  reduction  in  ohmic  losses  of  the  TEoi  mode  at  0.245  THz  when  R  grows  from  the  waveguide 
cut-off  value,  with  the  radii  of  interest  examined  in  the  main  figure  marked  as  square  markers. 
For  small  radius,  the  ohmic  loss  is  high  and  hence,  possesses  great  influence  on  the  direction 
of  threshold  shift.  In  this  regime,  when  the  radius  increases,  the  corresponding  fall  in  ohmic 
loss  has  the  tendency  to  lower  the  threshold.  On  the  other  hand,  at  larger  R  values  where  the 
ohmic  losses  are  low,  the  vibrational  bottleneck  increasingly  raises  or  prevents  threshold  with 
increasing  R  and  pressure.  However,  at  pressures  less  than  100  mTorr  in  large  cells  (i.e.  R  >  0.26 
cm),  the  pressure-dependent  threshold  lines  in  Fig.  4(c)  converge  to  the  linear  behavior  observed 
for  cells  with  R  <  0.26  cm.  As  we  shall  see  next,  this  follows  directly  from  the  fact  that  Amfp 
grows  large  at  low  pressures,  so  for  a  constant  intensity  and  pump  rate  Rpump,  the  inversion 
(AN)  and  gain  exhibit  little  radial  dependence.  As  a  result,  variation  of  R  has  a  much  smaller 
effect  on  the  wall  collision  rate,  and  correspondingly  on  (7),  at  lower  pressures.  So,  for  cell  radii 
where  the  ohmic  losses  are  also  low,  the  threshold  remains  almost  unchanged  as  R  is  varied. 

To  understand  the  spatial  effects  on  laser  performance  more  clearly,  the  radial  dependence 
of  the  inversion  and  gain  in  a  R  =  0.5  cm  cavity  pumped  by  CW  CO2  laser  at  Rpump  =  100  W 
are  studied  in  Fig.  5.  From  the  intensity-weighted  averaging  scheme  used  in  Eq.  (4),  we  have 
already  noted  that  it  is  most  desirable  to  match  the  peak  of  the  cavity’s  eigenmode  intensity 
to  regions  where  AN(r)  or  7(r)  is  maximum.  Such  consideration  plays  an  important  role  in 
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Pressure,  mTorr 


Fig.  5.  (a)  Radial  spatial  variation  of  the  excited  inversion  in  V3  state  of  *^CH3F  for  a  range 
of  operating  pressure  from  50  to  450  mTorr.  0  cm  marks  the  center  of  the  cylindrical  cavity 
while  0.5  cm  marks  the  copper  wall.  All  four  panels  in  this  figure  assumed  CW  pumping  at 
Fpump  =  100  W.  Near  the  cell  wall,  the  diffusion  is  ballistic  allowing  significant  reduction 
of  molecules  from  the  V3  pool,  and  hence,  results  in  a  higher  population  inversion,  (b)  Same 
as  in  (a)  except  that  the  transition  width  at  0.245  THz  is  also  factored  in  to  study  the  gain, 
(c)  Gain  dependence  on  pressure  for  a  fixed  set  of  pump  parameters.  From  a  molecular  gas 
physics  standpoint,  results  clearly  indicate  that  small  R  cavities  are  favored  in  terms  of  gain 
magnitude  and  pressure  cut-off.  (d)  Data  extracted  from  (c)  where  the  left  axis  and  right 
axis  respectively  depict  the  optimum  operating  pressure  and  the  corresponding  gain,  as  a 
function  of  cell  radius  R.  Again,  results  here  favored  small  sized  cavity  for  high  pressure 
operation. 


designing  the  cavity  itself.  Figure  5(a)  illustrates  the  distribution  of  AN{r)  within  the  cavity 
for  values  of  pressure  ranging  from  50  mTorr  to  450  mTorr.  In  general,  strong  inversions  are 
favored  near  the  cell  wall  where  the  wall  collision  rates  are  the  highest,  and  so,  allow  rapid 
depopulation  of  the  excited  and  higher  excited  thermal  pools  to  reduce  pool  quenching.  This 
directly  implies  an  advantage  in  exciting  a  lasing  mode  whose  intensity  peaks  near  the  cell 
wall.  However,  such  approach  will  also  have  to  be  weighted  against  the  higher  ohmic  losses 
that  arise  from  the  field’s  enhanced  interaction  with  the  cell,  in  a  similar  trade-off  effects  to 
what  was  already  discussed  in  Fig.  4(c).  The  details  and  treatment  of  the  cavity  design  will  be 
left  for  future  work.  We  further  note  that  the  magnitude  of  radial  variation  in  the  cell  increases 
with  pressure.  For  instance,  at  50  mTorr,  a  good  proportion  of  the  cavity  is  ballistic  [see  Amfp 
variation  in  Fig.  4(c)]  so  that  AN  is  approximately  uniform,  i.e.  AN (rcen)  =  0.96AA (rwaii)  where 
Teen  and  r^aii  indicate  the  position  at  the  cell  center  and  wall,  respectively.  On  the  other  hand,  at 
450  mTorr,  most  of  the  molecules  are  in  the  non-ballistic  regime  where  k^^\\  oc  1  jr^  as  opposed 
to  kuirect  1  /r.  This  leads  to  AA(rcen)  to  be  only  14%  of  AAt(rwaii)-  Indeed,  it  was  found  that 
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the  maximum  inversion  for  the  set  of  parameters  used  in  this  example  increases  with  pressure 
because  of  the  corresponding  increase  in  A^tot-  Figure  5(b)  depicts  a  similar  plot  for  the  gain. 
The  same  characteristics  are  observed  among  the  two  except  that  7(r)  is  also  influenced  by 
the  laser  transition  linewidth  which  broadens  as  pressure  increases.  This  decrease  of  gain  with 
pressure  is  again  consistent  with  the  discussion  presented  for  Fig.  4(c),  explaining  why  most 
OPFIR  lasers  use  large  diameter  cavities  operating  at  low  pressure. 

However,  this  design  misses  a  critical  point.  Since  Rpump  increases  when  the  same  pump 
power  is  more  tightly  confined,  the  gain  tends  to  grow  as  cavity  radii  shrink  until  it  eventually 
saturates  (because  A^tot  is  a  conserved  quantity).  Figures  5(c)  and  5(d)  track  the  optimum  radius 
size  for  lasing  based  on  the  numerical  model  developed.  An  optimum  operating  pressure  exists 
for  each  radius  where  the  gain  is  maximized.  Beyond  that,  vibrational  bottleneck  begins  to  play 
a  dominant  role,  explaining  the  fall  of  the  gain  as  pressure  rises.  At  high  enough  pressure,  the 
lasing  action  will  be  quenched.  For  example,  in  Fig.  5(c)  the  cut-off  pressure  occurs  at  380 
mTorr  for  R  =  I  cm.  Because  of  the  role  of  higher  lying  vibrational  levels,  the  cutoff  pressure 
is  pump  intensity  dependent  such  that  the  greater  the  pumping  intensity,  the  higher  the  cut-off 
pressure.  Thus,  the  optimum  pressure  and  gain  results  summarized  in  Fig.  5(d)  based  on  the 
molecular  gas  physics  greatly  favor  small-sized  cavities,  especially  if  the  effects  of  waveguide 
ohmic  losses  can  be  overcome. 

4.  Conclusions 

We  have  presented  a  comprehensive  theoretical  analysis  of  lasing  action  in  optically-pumped 
rotationally  excited  molecular  gases  confined  in  a  metallic  cavity.  In  order  to  do  that,  we  have 
developed  a  theoretical  framework  that  provides  a  realistic  spatio-temporal  description  of  the 
molecular  de-excitation  and  diffusion  processes  occurring  inside  the  cavity.  Our  approach  also 
accounts  for  the  influence  on  the  main  laser  features  of  the  electric-field  intensity  distribu¬ 
tion  and  the  ohmic  losses  associated  to  the  different  cavity  modes.  To  illustrate  this  theoretical 
framework,  we  have  analyzed  in  detail  the  case  of  lasing  action  from  optically -pumped  methyl 
fluoride  gas  confined  in  a  cylindrical  copper  cavity.  Our  simulations  have  shown  that  for  this 
system,  the  region  of  maximum  population  inversion  is  located  near  the  cavity  walls.  The  out¬ 
come  of  our  calculations  has  also  shown  that  the  lowest  lasing  threshold  mode  corresponds  to 
the  cavity  mode  whose  electric-field  intensity  profile  presents  an  optimal  spatial  overlap  with 
the  corresponding  population  inversion  distribution  inside  the  cavity,  and  at  the  same  time,  mi¬ 
nimizes  the  ohmic  losses  at  the  lasing  frequency.  Furthermore,  we  have  predicted  the  existence 
of  an  optimal  radius,  determined  by  the  interplay  between  the  overall  gain  in  the  system  and 
the  corresponding  ohmic  losses,  which  minimizes  the  lasing  threshold  intensity  for  a  large  in¬ 
terval  of  gas  pressures.  We  believe  the  results  reported  in  this  paper  could  contribute  to  the 
development  of  a  novel  class  of  room-temperature  compact  sources  of  far-infrared  and  THz 
radiation. 

Appendix:  Rate  equations  and  rate  constants  of  ^^CH3F  gas  lasers 

In  this  appendix,  we  first  present  the  set  of  rate  equations  used  to  model  the  *^CH3F  gas  lasers. 
Figure  6  and  Eq.  (7)  to  Eq.  (15)  refer  to  a  specific  example  of  a  methyl  fluoride  gas  discussed 
in  Sec.  3,  adopted  from  the  general  model  presented  in  Eig.  1,  Eq.  (1)  and  Eq.  (2)  of  Sec.  2. 
We  have  labeled  the  rotational  levels  from  Ni  to  Ag,  and  the  thermal  pools  as  At,  Nb  and 
Ac.  Hence,  Tn,  T23,  T45,  T56,  and  their  reverse  processes  represent  rotational  transitions  due  to 
dipole-dipole  collisions,  T\a,  T4b,  and  similar  processes  represent  K-swap  transition  lifetimes 
from  the  rotational  levels  (i.e.  Ai  and  A4)  into  their  respective  thermal  pools  (i.e.  Na  and  Nb), 
and  Tab,  '^ac,  and  their  reverse  processes  are  vibrational  transitions  which  arise  due  to  collisions 
with  the  cell  wall.  Readers  should  note  that  all  notations  in  the  main  body  of  the  paper  are  kept 
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Fig.  6.  *^CF[3F  OPFIR  laser  model  used  in  this  paper:  results  of  diagnostic  studies  and 
theoretical  considerations  restrict  the  degrees  of  freedom  so  that  a  numerically  tractable, 
yet  physically  accurate  model  is  attained.  Details  on  each  process  labeled  are  described  in 
the  text. 


consistent  with  those  used  in  Fig.  1,  Eq.  (1)  and  Eq.  (2). 
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where  AA^pump(r,f)  =  Nf''(r,t)  —  g5N^'-{r,t)/g2  with  N^'-(r,t)  being  the  total  molecular  den¬ 
sity  of  7  =  4  nonthermal  level  in  Vq  state  after  accounting  for  the  fractional  contribution  from 
thermal  pool  A,  and  g2  =  9  is  its  degeneracy.  Similar  definitions  hold  for  and  gs  =  1 1 . 

Table  1  provides  the  rate  constants  and  relevant  molecular  parameters  used  in  our  calculations 
[8,14,15], 


Table  1.  Translational  Cross  Sections  (in  of  *^CH3F 


Mass  (AMU) 

35 

Gas  kinetic  collision  cross  section 

44 

V3  —1"  W  cross  section  at  300  K,  ffexc 

2.65 

Dipole-dipole  cross  section  at  300  K,  Odd 

320 

K-swap  cross  section  at  300  K,  Ok-swap 

137 
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